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Introduction 


This  study  mostly  considers  three  basic  classes  of  free-surface  flows  related  to  interactions  of 
water  waves  with  floating  bodies  like  ships  and  offshore  structures  :  (i)  time-harmonic  flow  associ¬ 
ated  with  diffraction-radiation  of  regular  waves  by  a  floating  body  without  forward  speed,  e.g.  an 
offshore  structure,  (ii)  steady  flow  due  to  a  ship  advancing  in  calm  water,  (iii)  time-harmonic  flow 
due  to  a  ship  advancing  in  regular  waves,  i.e.  wave  diffraction-radiation  with  forward  speed.  The 
body  of  water  is  assumed  to  be  of  infinite  depth  and  lateral  extent. 

The  classical  potential-flow  theoretical  framework  is  adopted.  This  framework  correctly  repre¬ 
sents  the  aspects  of  ship  wavemaking  and  wave  diffraction-radiation  by  ships  and  offshore  struc¬ 
tures  that  are  of  primary  practical  significance  (notably  forces,  pressure  distributions,  and  wave 
patterns) .  Furthermore,  potential-flow  theory  is  the  only  framework  that  provides  practical  means 
of  analyzing  steady  ship  waves  and  wave  diffraction-radiation  in  various  wave  environments  for 
routine  applications,  and  for  hull-form  optimization  (ultimately  required  for  superior  designs) . 

Mathematics  has  been  used  extensively  to  analyze  free-surface  flows  (and  more  generally  fluid 
flows)  in  various  special  cases,  including  special  geometries  (e.g.  2D  flows,  and  3D  flows  about  simple 
bodies  like  spheres  and  spheroids)  and  limiting  cases  associated  with  a  perturbation  analysis  (e.g. 
thin-ship  and  slender-body  theories,  and  theories  based  on  low  or  high  ship-speed  or  wave-frequency 
approximations) .  This  traditional  role  of  mathematics,  which  may  be  named  "mathematics  for 
special  cases”,  has  provided  a  wealth  of  useful  results.  Another  important  role  of  mathematics, 
which  may  be  named  "mathematics  for  the  computer  age”,  consists  in  developing  mathematical 
formulations  that  make  it  possible  to  perform  robust,  accurate,  and  efficient  numerical  calculations, 
without  restriction  with  regard  to  hull  form,  ship  speed,  and  wave  frequency.  Mathematics  is  used 
here  for  the  latter  purpose. 

Specifically,  four  fundamental  analytical  representations  of  the  three  classes  of  free-surface  flows 
under  consideration  are  summarized.  These  analytical  representations  are  the  main  results  that 
underly  the  Fourier-Kochin  theory  of  ship  waves.  The  four  analytical  representations  are  presented 
in  an  order  and  a  manner  intended  to  provide  a  coherent  and  complete  overview  of  this  theory. 
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Two  complementary  elementary  solutions  of  Laplace’s  equation 
An  important  fundamental  solution  of  the  Laplace  equation  is  the  elementary  wave  function 

exp[  \J  o?-  +0"^  C  +  0'n)]  (la) 

Here,  a  and  (i  are  real.  Another  important  class  of  fundamental  solutions  of  Laplace’s  equation  is 
the  elementary  Rankine  singularity  \/r  and  the  related  Rankine  sources  1/r'  and  1/r"  with  r,  r', 
r"  given  by 

r  =  ^/{i-xY+{r]-yY+{(i-zY 

r'  =  ^/  {r)-yY+  (C+'2)'^  ^  (^'^) 


Here,  (7^  is  a  positive  real  constant. 

The  elementary  wave  (la)  and  the  elementary  Rankine  singularity  1/r  are  related  via  the  Fourier 
transformation  (see  Appendix  1) 

S 


1  I  too  roo 


(2a) 


with 


k  =  £:  =  e-i[a(C-x)  +  /?(r,-y)]  (2b) 

Expression  (2a)  shows  that,  for  ^  +  2  <  0 ,  the  related  Rankine  sources  1/r'  and  1/r"  can  similarly 
be  expressed  as 

The  Fourier  relations  (2)  indicate  that  free-surface  flows  can  in  principle  be  represented  using  ex¬ 
clusively  Rankine  singularities  or  elementary  waves.  However,  these  two  "extreme”  approaches  are 
not  the  most  appropriate  because  Rankine  singularities  and  elementary  waves  are  complementary 
elementary  solutions  of  the  Laplace  equation  that  are  ideally  suited  for  representing  nonoscillatory 
local  flow  disturbances  and  waves,  respectively. 
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Free-surface  boundary  conditions 


Consider  a  ship  advancing  in  waves  along  a  straight  path,  with  constant  speed  U .  The  X  axis 
is  chosen  along  the  path  of  the  ship,  and  points  toward  the  ship  bow.  Thus,  the  ship  advances  in 
the  direction  of  the  positive  X  axis.  The  Z  axis  is  vertical  and  points  upward,  and  the  mean  free 
surface  is  taken  as  the  plane  Z  =  0 .  The  flow  is  observed  from  a  Cartesian  system  of  coordinates 
moving  with  speed  U  along  the  X  axis,  and  is  expressed  as  the  sum  of  a  steady-flow  component, 
which  represents  the  flow  due  to  the  ship  advancing  in  calm  water,  and  an  unsteady-flow  component 
associated  with  the  ambient  waves,  the  waves  diffracted  by  the  ship,  and  the  radiated  waves  due 
to  the  motions  of  the  ship  about  a  mean  position. 

Within  the  framework  of  a  linear  potential-flow  analysis,  the  velocity  potential  $  associated 
with  the  unsteady  flow  due  to  the  ship  satisfies  the  free-surface  boundary  condition 

g^z  +  {dr  -  Udx)'^  ^  =  0  at  Z  =  0 

where  T  stands  for  time  and  g  is  the  gravitational  acceleration.  A  frequency-domain  analysis  is 
considered.  Thus,  the  time-harmonic  flow  due  to  a  ship  advancing  through  an  elementary  wave  is 
now  examined.  The  velocity  potential  of  a  linear  time-harmonic  flow  can  be  defined  as  the  limit, 
when  the  growth  parameter  ^  ►  -|-0 ,  of  the  flow  given  by  the  potential 

3ie  $(X)  exp(/xr/2-iwr) 

where  3ie  stands  for  the  real  part  and  w  is  the  frequency  of  the  waves  encountered  by  the  advancing 
ship.  This  expression  for  the  potential  is  associated  with  a  flow  that  may  be  imagined  to  start  at 
the  time  T  —  — oo  from  a  state  of  rest  (since  /z  >  0) ,  corresponding  to  the  initial  conditions  $  =  0 
and  =  0 .  The  free-surface  boundary  condition  becomes 

g^z  —  +*m/2  —  iUdx)'^^  =  ^ 

The  time-harmonic  flow  due  to  the  ship  can  be  formulated  in  nondimensional  form  in  terms  of  the 
water  density,  a  reference  length  L  and  a  reference  velocity  U .  E.g.,  L  and  U  may  respectively  be 
taken  as  the  length  and  the  speed  U  of  the  ship.  The  free-surface  condition  associated  with  the 
nondimensional  coordinates  ^  =  X/L  and  potential  (j)  =  ^/{UL)  is 

<j)^  —  (f +  ie/2  —  i  F  8^)“^  ^ 

with 

f  =  ujyfL/g  e^nyflfg  F  =  Ul^f^ 
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For  0  <  e  <C  1 ,  the  free-surface  condition  becomes 


4>c,-{f  4)-ie{f  -iFd^)4>  =  {) 

Thus,  the  linearized  free-surface  boundary  condition  associated  with  diffraction-radiation  of 
time-harmonic  waves  with  forward  speed  is 

(t>Q  -  f<t>  +  i 2r  <t>^  -f  -£{if(l>  +  F<f)^ )  =  p  (3a) 

where  t  =  fF  is  the  Brard  number,  and  p  stands  for  a  nonhomogeneous  forcing  term.  E.g.,  p  may 
account  for  a  pressure  distribution  at  the  free  surface.  The  linear  free-surface  condition  associated 
with  steady  ship  waves  can  similarly  be  shown  to  be 

</><  +  c  -  eF(p^  =  p  (3b) 

Here,  stands  for  the  disturbance  velocity  due  to  the  ship.  Finally,  the  free-surface  condition  for 
linear  time-harmonic  flows  without  forward  speed  is 

<f>C  -  f(t>-i£N  =  P  (3c) 

The  free-surface  conditions  (3b)  and  (3c)  evidently  correspond  to  the  special  cases  /  =  0  and  F  =  0 
of  the  ship)-motion  condition  (3a) . 


Dispersion  functions  and  dispersion  curves 


The  elementary  wave  solution  (la)  of  the  Laplace  equation  satisfies  the  homogeneous  free-surface 
boundary  conditions  (3)  if  the  Fourier  variables  a  and  0  satisfy  the  relation  Dg  =  0,  where  the 
function  =  0  is  defined  as 

Dg  =  D  i  £  D'  (4a) 


with  D  and  D'  given  by 


'{f-Fa)^-k' 

f-Fa 

/F>o' 

F'^a^  -  k 

>  D'  =  < 

—Fa 

>  if  ^ 

f  =  0  > 

f-k 

\  > 

/ 

k.  s 

(4b) 


In  expressions  (4b)  for  the  function  D,  called  dispersion  function,  k  =  is  the  wavenum¬ 

ber.  Farfield  waves  corresponding  to  the  free-surface  boundary  conditions  (3)  are  given  by  one¬ 
dimensional  Fourier  superpositions  of  elementary  waves  (la)  since  the  Fourier  variables  a  and  0 


4 


are  related  by  the  dispersion  relation  D  =  0.  In  fact,  the  dispersion  relation  £>  =  0  typically  defines 
several  curves  (although  a  single  curve  may  be  defined  in  some  simple  cases)  in  the  Fourier  plane 
{a,0) .  These  curves,  called  dispersion  curves,  have  an  important  role  in  the  Fourier  representation 
of  farfield  and  nearfield  waves.  The  dispersion  curves  associated  with  the  dispersion  relations  (4b) 
are  now  considered. 

For  time-harmonic  flows  without  forward  speed,  the  dispersion  relation  £>=0  defines  the  single 
closed  dispersion  curve 

k  =  f  (5a) 

which  evidently  is  a  circle  centered  at  the  origin  A:  =  0  of  the  Fourier  plane.  For  steady  flows,  the 

dispersion  relation  k  =  defines  two  open  dispersion  curves  given  by 

(3  =  ±\a\  (5b) 

with  F^|a|  >  1 .  These  dispersion  curves  are  symmetric  about  both  f3  =  0  and  a  =  0 .  The  ship- 
motion  dispersion  relation  D  =  0  corresponding  to  the  case  r  =  fF  >  0  yields 

/?  =  ±^{f-FaY-a^  (5c) 

The  dispersion  curves  defined  by  (5c)  are  symmetric  with  respect  to  /3  =  0.  The  ship-motion 
dispersion  relation  (5c)  defines  two  or  three  dispersion  curves  if  r  is  larger  or  smaller  than  1/4 . 

If  T  <  1/4 ,  the  three  dispersion  curves  intersect  the  axis  P  =  0  &t  four  values  of  a ,  which  are 
denoted  af  and  and  are  given  by 

'±af  =  /2(V!74±T-1/2)Vt^' 

<  _  >  (6) 

±a±  =  (v/l/4±r  +  ll2flF^ 

V  -- 

The  three  dispersion  curves  are  found  in  the  regions 

a  <a~  OL~  <a  <Q^  Q: <  o;  (7a) 


If  T  =  1/4 ,  (6)  yield 

a-=aT  =  -f/F=-uLIU 

and  the  dispersion  curves  in  the  regions  a <  and  a~  <ol<  oil  connected.  If  r  >  1/4 ,  the 
roots  a~  and  a~  are  complex  and  we  only  have  two  dispersion  curves,  located  in  the  regions 

a  <al  al<a  (7b) 
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The  inner  roots  af  and  the  outer  roots  af  satisfy  the  inequalities 

{a;<-f/F<a-<0<at<f/F<at] 

The  dispersion  curves  defined  by  (5),  (6)  and  (7)  are  depicted  for  r  =  0.2 , 0.25 , 0.3  in  Fig.  1 ,  where 
the  Fourier  plane  is  scaled  with  respect  to  f  /F . 

In  the  limit  r  =  0,  (6)  yield  —  and  a*  =  ±1/F^.  More  generally,  we  have 

w  ±/^  af^±l/F^  forr<l 

Thus,  the  dispersion  curves  in  the  outer  regions  a  <  a~  and  aj  <  a  correspond  to  values  of 
the  wavenumber  k  that  are  greater  than  1/F^  (approximately),  whereas  we  have  p  for  the 
dispersion  curve  in  the  inner  region  Ci~  .  The  ratio  of  the  wavenumbers  ki  and  ko  associated 

with  the  inner  and  outer  dispersion  curves,  given  by  ki/ko  ~  ,  is  very  small  if  r  is  small,  say  if 

t<0.2  .  Thus,  scaling  of  the  Fourier  coordinates  a  and  jS  with  respect  to  p  and  1/F^  is  proper  in 
this  regime.  The  upper  and  lower  parts  of  Fig.  2  respectively  depict  the  inner  and  outer  dispersion 
curves  for  several  values  of  r  in  the  frequency-scaled  and  speed-scaled  Fourier  planes  {a  ,(i)/ P  and 
F^{a,P).  Finally,  (6)  yield 

af^f/F^Og  asr— >'00 

Thus,  scaling  of  the  Fourier  coordinates  with  respect  to  //F  is  proper  for  large  values  of  r . 

Free-surface  Green  functions 

The  Green  function  G{^,x)  associated  with  the  ship-motion  free-surface  boundary  condition 
(3a)  satisfies  the  free-surface  boundary  condition 

Gc  -pG  +  i  2tG^  +  F^G^^  -e{ifG  +  FG^)  =  0  (9) 

atC  =  0.  G{^,x)  is  the  velocity  potential  of  the  flow  created  at  the  field  point  ^  =  (^ , ?/ ,  C ^ 0) 
by  a  pulsating  source  located  at  the  singular  point  x  =  (x  ,y  ,z  <0) .  Green  functions  associated 
with  a  boundary  condition  like  (9) ,  which  holds  at  the  plane  C  =  0  >  can  be  expressed  as 

G({,  X )  =  -x,Tf-y,C-z)+  G^(^ -x,Tj-y,C  +  z)  (10a) 

where  G^  corresponds  to  a  simple  Rankine  source  and  G^  accounts  for  free-surface  effects.  The 
Rankine  component  G^  is  defined  by 

47rG^=— 1/r  (10b) 
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with  r  given  by  (lb) .  The  free-surface  component  in  expression  (10a)  corresponding  to  the 
free-surface  condition  (9)  is  given  by  (see  Appendix  2) 

1 

47rG'^=  lim  -  d0  da^  +  (11a) 

e^+o  TT  J-oo  J-00  De  ^  ’ 

Here,  k  and  £  are  given  by  (2b),  Dg  is  the  dispersion  function  defined  by  (4) ,  and  the  amplitude 
function  is  given  by 

=  [{f-Fa)^+k]/{2k)  (11b) 


The  free-surface  component  in  (10a)  is  usually  expressed  in  the  alternative  forms 


47rG^  = 


-l/r'+H 
<  ► 

l/r'+H 


(12) 


with  r'  given  by  (lb).  The  alternative  decompositions  47rG'^=  —Ijr'+T-l  and  47r(j^=  \/r'+H 
are  shown  further  on  to  be  appropriate  for  F=0  and  F>0,  respectively.  Thus,  (10)  and  (12)  yield 


47rG  =  ^ 

-l/r  -  Ijr'+H 

1  1 

>  if  < 

V  =  0 

> 

\-l/r  +  l/r'+H\ 

[f>o| 

(13a) 


Expressions  (2c)  and  (11)  show  that  the  free-surface  components  Ti  and  H  in  (13a)  are  given  by 
the  Fourier  superposition  of  elementary  waves 


1  roo  roo  A 

jr=  lim  -  dl3  + 

£~^+0  TT  J —oo  J —oo 


(13b) 


with  Dg  given  by  (4)  and 


F=l 


n 

H 


A  = 


Plk\ 


(13c) 


In  the  special  case  of  wave  diffraction-radiation  without  forward  speed,  (13c)  yields  A  =  liik  =  p, 
i.e.  at  the  dispersion  curve  defined  by  the  dispersion  relation  D  —  Q. 
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Two  fundamental  boundary- integral  representations 

Consider  the  basic  problem  of  determining  the  flow,  inside  a  free-surface  potential-flow  region, 
that  corresponds  to  a  given  flow  at  a  surface  E  bounding  the  potential-flow  region.  Here,  the 
boundary  surface  S  stands  for  E  =  E^U  where  the  surface  E^  is  an  arbitrary  control  surface 
outside  the  viscous  boundary  layer  that  surrounds  the  hull  of  a  ship  (or  other  body,  e.g.  an  offshore 
structure)  at  or  below  the  free  surface,  and  the  surface  E^  represents  the  outer  edge  of  the  viscous 
wake  trailing  the  ship  or  a  control  surface  outside  the  viscous  wake.  If  viscous  effects  are  ignored, 
E^  may  be  taken  as  the  wetted  surface  of  the  ship  hull.  For  a  ship  equipped  with  lifting  surfaces, 
e.g.  a  sailboat,  E*^  includes  the  two  sides  of  every  vortex  sheet  behind  the  ship.  For  a  multihull 
ship,  the  hull-l-wake  surface  E  consists  of  several  component  surfaces,  which  correspond  to  the 
separate  hull  components  of  the  ship  and  their  wakes. 

E  may  intersect  the  mean  free-surface  plane  ^  =  0  or  may  be  a  closed  surface  entirely  located 
below  the  plane  C  =  0  •  In  the  former  case,  E  is  bounded  upward  by  its  intersection  curve  F  with 
the  plane  C  =  0 .  The  unit  vector  n=  ( n^,  n^,  )  normal  to  the  boundary  surface  E  points  outside 

E ,  i.e.  inside  the  flow  domain.  The  unit  vector  t  =  (t^,t^,0)  tangent  to  the  boundary  curve  F  is 
oriented  clockwise  (looking  down).  The  unit  vector  normal  to  F  in  the  plane  C  =  0  points  inside 
the  flow  domain,  like  ft ,  and  is  given  by  i7  =  (-<*',  t®,  0 ) . 

A  classical  fundamental  result  of  potential-flow  theory  is  that  the  velocity  potential  (f)  within 
a  potential-flow  region  bounded  by  a  surface  E  is  explicitly  defined  in  terms  of  the  values  of  the 
potential  4>  and  its  normal  derivative  d4>/dn  =  V<f)  ■  n  =  u  •  n  &t  the  boundary  surface  E .  This 
fundamental  boundary-integral  representation  —  called  potential  representation  because  it  defines 
the  potential  (f)  —  and  related  boundary-integral  representations  that  define  </>  within  a  potential- 
flow  region  in  terms  of  boundary  distributions  of  sources  or  normal  dipoles  form  the  basis  of 
numerous  panel  methods  reported  in  the  ship  and  offshore-structure  hydrodynamics  literature  to 
compute  linear  free-surface  potential  flows  using  free-surface  Green  functions. 

Another  fundamental  boundary-integral  representation  —  given  in  Noblesse  (2000)  and  called 
velocity  representation  —  explicitly  defines  the  velocity  u  within  a  potential-flow  region  in  terms 
of  the  velocity  u  at  the  boundary  surface  E .  Thus,  the  velocity  representation  does  not  involve 
the  potential  cf> .  This  property  is  a  major  difference  between  the  velocity  representation  and  the 
potential  representation.  The  velocity  representation  can  therefore  be  used  to  couple  a  viscous 
flow,  for  which  a  velocity  potential  cannot  be  defined,  and  a  potential  flow  in  a  direct  manner,  i.e. 
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without  having  to  solve  an  integral  equation  as  is  necessary  if  the  potential  representation  is  used. 
Specifically,  the  potential  representation  requires  a  two-step  procedure  :  an  integral  equation  must 
first  be  solved  to  determine  <j)  at  the  boundary  surface  from  the  known  boundary  value  of  d(l)/dn , 
and  0  and  can  subsequently  be  determined  at  interior  points  from  the  boundary  values  of  (f) 
and  d(t)/dn. 

Another  important  difference  between  the  potential  representation  and  the  velocity  representa¬ 
tion  is  that  the  velocity  representation  defines  u  in  terms  of  first  derivatives  of  the  Green  function 
G ,  whereas  u  can  only  be  obtained  from  the  potential  representation  using  either  analytical  differ¬ 
entiation  of  the  potential  representation,  which  involves  second-order  derivatives  of  G ,  or  numerical 
differentiation  of  (j). 

The  velocity  representation  and  the  classical  potential  representation  are  given  in  Noblesse 
(2000)  for  generic  free-surface  potential  flows  and  for  four  basic  classes  of  flows  corresponding  to  (i) 
an  infinitely  rigid  or  soft  free-surface  plane,  (ii)  diffraction-radiation  of  time-harmonic  water  waves 
without  forward  speed,  (iii)  steady  ship  waves,  and  (iv)  time-harmonic  ship  waves  (diffraction- 
radiation  with  forward  speed).  The  classical  potential  representation  and  the  velocity  representation 
are  given  here  for  diffraction-radiation  of  time-harmonic  water  waves  without  forward  speed,  steady 
ship  waves,  and  time-harmonic  ship  waves.  The  potential  4>  and  the  velocity  u  can  be  expressed  as 


(j)  =  (f)^  4-  (j)^  +  (fp^  —  (j)P 
-t-  -I- 


^  if  F  =  0 


^  =  (j)^  -  (j)^  +  (f)^  +  (j)^  - 

u  =  -f  U^+  Up  — 


S  if/  =  0 


(14a) 


(14b) 


<i)  -  (j)R  +  +  (t>f  +  (jiH  -  (i>P 

U  =  +  U^+  Up  +  u^  +  u^  —  u^ 


if  /F  >  0 


(14c) 


The  components  0^,  uR  and  0;^,  u^  in  the  potential  and  velocity  representations  (14)  stem 
from  the  Rankine  components  in  the  Rankine-Fourier  decompositions  (13a)  of  the  Green  functions 
associated  with  the  free-surface  boundary  conditions  (3) .  The  Rankine  components  in  (14)  are 
given  by  surface  distributions  of  elementary  Rankine  singularities  1/r  and  1/r'.  Specifically,  the 
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Rankine  component  (p^,  is  given  by 


'-(pf^ 

—{u-n)R  +  (p  {n^Rx  +  n^Ry  +  n^Rz) 

>=  d A  < 

-{u-  n)  Rx  —  {uxfiy  Ry  +  (uxn)*'  Rz 

yR 

Jy. 

—{u-n)Ry  —  {uxfiy  Rz  +  {uxfiy  Rx 

—{u‘fi)Rz  -  {uxfiy Rx  +  {uxfiy  Ry  ^ 

with  R  —  —If r .  The  Rankine  component  is  given  by 


'-<P^' 

' -{u-fi)R* +  (P  (n^R*  +  nyRl  +  n^Rl) ' 

47r  < 

>  =  dA< 

-(u-n)R*  -  {uxfiy  R*  +  {uxfiy  Rl 

JT, 

—{u-fi)Ry  —  {uxfiy  Ri  +  {uxfiy  Ri 

—w^ 

-{u-fi)R*  -  {uxfiy  Rl  +  {uxfiy  R* 

(15a) 


(15b) 


with  R*  =  -l/r' . 

The  components  (p'^,  and  (p^ ,  in  the  potential  and  velocity  representations  (14)  are 

defined  in  terms  of  the  functions  H  and  H ,  which  account  for  free-surface  effects.  The  free-surface 
component  ^  —  where  T  stands  for  either  H  or  H  —  is  given  by 


'-<p^' 

-{u  •n)T-¥(p  {n^Tx  +  vi}>Ty  +  n^Tz) 

u^ 

>  =  dA< 

—{u  •fi)!Fx  —  {uxfiy  Ty  +  {uxfiy  Tz 

Je 

—{u  •n)Ty  —  {uxfiy  Tz  +  {uxfiy  Rx 

—w^ 

—{u  ■n)Rz-  {uxfiy Rx  +  {uxfiy Ry 

Expressions  (15)  and  (16)  involve  source  and  vortex  distributions  over  the  boundary  surface  E 
with  respective  strengths  given  by  the  normal  component  u-n  and  the  tangential  component  uxn 
of  the  velocity  at  E.  The  velocity  components  in  (15)  and  (16)  only  involve  the  first  derivatives 
VR,  VR*,  VR. 

The  component  in  (14a) ,  the  component  (pp  ,  Up  in  (14b)  and  (14c) ,  and  the  components 
(pj  ,  uf  and  uj^  in  (14c)  are  given  by  distributions  over  the  boundary  curve  T .  Specifically,  the 
component  uj^  in  (14a)  is  given  by 


f  > 

u^ 

t^Rz 

Air  < 

vj^ 

>  =  j  dC  (p  < 

-e  Hz  > 

wy 

K  J  J 

tyHx-t^Hy 

K  ^  / 

(17a) 
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The  component  ,  Up  in  (14b)  and  (14c)  is  given  by 


'  {u-u)  {ty)^H  -{u-t)  etyR  +  tyR^^ ' 

^  =  F^j^dC  < 

{u-P){ty)^R^-{u-t)tHyR:c 

(u-u)  {ty)^Ry  -iu-t)  (tnyRy+R^) 

(u-p)  {ty)^R,  -(u-t)  {tnyR.-R^^y)  ^ 

The  components  <j)^ ,  and  in  (14c)  are  given  by 


(t)tyR 

u? 

f 

0 

< 

>  =  2  2t  /  dC< 

Jr 

u-t  R 

w? 

u-t  R^ 

/  N 

Uf 

tyR 

47r  ^ 

vf 

II 

to 

fdC(l)< 

r 

-t^R  ^ 

wH 

\  J  J 

tyR^-t^R^ 

\  y  J 

The  function  in  (17b)-(17d)  is  defined  as 


(17b) 


(17c) 


(17d) 


H^^-x,ri-y,C  +  z)  =  [  dt  H(^  -  x  ,t]  -  y +  t) 

The  function  H^y  in  (17b)  is  comparable  to  the  first  derivative  VH  of  H  (the  functions  H^y  and 
VH  have  similar  behaviors  in  the  near  field  and  the  far  field) .  The  functions  and  in  (17c) 
and  (17d)  likewise  are  comparable  to  the  function  H . 


Finally,  the  component  (j)^,  associated  with  the  nonhomogeneous  term  p  in  the  free-surface 
conditions  (3)  is  given  by 


f  \  f  \ 


4>p 

R 

R 

uP 

f  , . 

-Rx 

f 

He 

>=  dAp< 

>=  dAp< 

vP 

-Ry 

Rv 

J 

He 
\  ^  / 

where  stands  for  the  mean  free  surface. 


(18) 


The  velocity  components  Up  and  defined  by  (17b)  and  (17c)  involve  the  normal  component 
u-u  and/or  the  tangential  component  u-t  of  the  velocity  u  at  the  boundary  curve  T .  Expressions 
(17a)  and  (17d)  for  the  components  and  uj^  involve  the  potential  (f>,  which  is  defined  in  terms 
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of  the  velocity  u  by  the  free-surface  conditions  (3c)  and  (3a) .  In  the  special  cases  /  =  0  and  F  =  0 , 
(14c)  can  be  verified  to  be  identical  to  (14b)  and  (14a)  as  expected. 

The  boundary-integral  representations  (14)-(18)  are  obtained  in  Noblesse  (2000)  using  elemen¬ 
tary  identities  in  vector  calculus.  Specifically,  the  divergence  theorem  is  used  in  the  3D  flow  region 
bounded  by  the  closed  surface  E  U  E°°  U  E^.  Here,  E°“  represents  a  large  outer  boundary  surface 
that  surrounds  the  previously-defined  hull-fwake  inner  boundary  surface  E ,  and  E^  stands  for  the 
portion  of  the  mean  free-surface  plane  z  —  0  that  is  bounded  by  the  intersection  curves  F  and  r°° 
between  the  plane  z  =  0  and  the  inner  and  outer  boundary  surfaces  E  and  E°°.  Stokes’  theorem  is 
also  used  in  the  2D  free-surface  region  E^  bounded  by  F  U  F°°  for  the  purpose  of  expressing  the 
contribution  of  the  free  surface  in  terms  of  a  line  integral  around  the  inner  boundary  curve 
F .  Thus,  the  line  integrals  (17)  and  the  corresponding  components  in  (14)  are  associated  with  the 
contribution  of  the  mean  free-surface  plane  E^  outside  the  boundary  surface  E .  The  flow  repre¬ 
sentations  (14)-(18),  given  here  for  deep  water,  can  be  extended  to  uniform  finite  water  depth,  and 
indeed  can  be  extended  to  a  broader  class  of  problems  involving  dispersive  waves  that  propagate 
in  the  presence  of  a  planar  boundary. 

The  classical  Green-function  method 

The  potential  and  velocity  representations  (14)-(18)  provide  a  mathematical  basis  for  computing 
steady  ship  waves  and  wave  diffraction-radiation  by  offshore  structures  and  ships  using  the  free- 
surface  Green  functions  (13) .  Free-surface  effects  in  the  boundary-integral  representations  (14)-(18) 
and  the  corresponding  free-surface  Green  functions  are  represented  by  the  components  7i  and  H 
given  by  the  double  Fourier  integral  (13b) . 

The  special  case  F  =  0  is  sufficiently  simple  that  the  double  Fourier  integral  (13b)  can  be  easily 
evaluated.  Specifically,  the  function  H  is  expressed  in  Noblesse  (1982)  as 

H  =  W  +  C  (19a) 

where  W  and  C  respectively  represent  a  system  of  circular  waves  and  a  local  flow  disturbance.  The 
wave  component  W  is  given  by  the  remarkably  simple  expression 

W  =  -27r/^  [Eo{f‘^h)  +iJo{f^h)]  exp[/^(C  +  2)]  (19b) 

obtained  in  Noblesse  (1982) .  Here,  h  =  \/(^  -  x)'^+  [r]  -  yY  ,  Jo  is  the  usual  Bessel  function  and 
Eq  is  the  Weber  function  defined  in  chapter  12  of  Abramowitz  &  Stegun  (1965) .  The  functions  Eq 
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and  Jo  are  infinitely  differentiable  and  can  be  evaluated  quite  efficiently.  The  local  ffow  component 
C  in  (19a)  can  be  evaluated  in  a  highly  efficient  manner  using  a  simple  analytical  approximation 
based  on  the  nearfield  and  farfield  analytical  approximations  given  in  Noblesse  (1982) ,  or  more 
accurately  using  table  interpolation  as  in  Ponizy  et  al.  (1994) . 

Simple  analytical  expressions  like  (19b)  are  not  available  if  F>0 ,  i.e.  for  the  free-surface  Green 
functions  associated  with  steady  and  time-harmonic  ship  waves,  because  forward-speed  effects  intro¬ 
duce  considerable  mathematical  difficulties.  More  precisely,  the  classical  Green-function  approach 
requires  practical  methods  for  performing  two  necessary  basic  tasks:  (i)  evaluation  of  the  free- 
surface  component  H  in  (13a)  and  its  gradient  Vff ,  i.e.  evaluation  of  the  double  Fourier  integral 
(13b) ,  and  (ii)  integration  of  H  and  VH  over  a  hull  panel  and  a  waterline  segment.  Both  of  these 
two  basic  tasks  are  difficult  because  the  functions  H  and  VH  are  singular  and  highly  oscillatory  in 
the  limit  r'  ^  0 .  Although  the  evaluation  and  the  subsequent  hull-panel  integration  of  H  and  VH 
have  been  considered  in  numerous  studies  in  the  ship-hydrodynamics  literature,  no  simple  practi¬ 
cal  method  of  performing  these  two  basic  tasks  appears  to  be  available  yet.  Indeed,  the  integral 
representations  of  the  Green  functions  associated  with  the  free-surface  boundary  conditions  for 
steady  and  time-harmonic  ship  waves  given,  for  instance,  in  Noblesse  (1981)  and  Chen  (1999)  are 
relatively  complicated  and  ill-suited  for  efficient  flow  calculations.  Free-surface  effects  associated 
with  the  Fourier  component  H  defined  by  the  double  Fourier  integral  (13b)  can  be  analyzed  more 
simply  and  effectively  using  the  Fourier-Kochin  approach  expounded  in  Noblesse  and  Yang  (1995) 
and  Noblesse  (2000) .  This  approach  is  summarized  below. 

Fourier-Kochin  representation  of  free-surface  eflfects 

The  order  in  which  the  tasks  of  evaluating  the  functions  H  and  VH  (a  Fourier  integration)  and  of 
integrating  H  and  VH  over  hull  panels  and  waterline  segments  (a  space  integration)  are  performed 
is  reversed  in  the  Fourier-Kochin  approach.  Specifically,  expression  (13b)  for  the  free-surface  com¬ 
ponent  H  is  substituted  into  (16)-(18) ,  and  the  integrations  with  respect  to  the  space  variables 
{x,y  ,z)  and  the  Fourier  variables  (a ,  /3)  are  permuted.  Thus,  space  integration  is  performed  first, 
and  Fourier  integration  is  performed  subsequently  in  the  Fourier-Kochin  approach. 

An  evident  advantage  of  the  Fourier-Kochin  approach  is  that  the  space  integration,  which  is  very 
complex  within  the  classical  Green-function  method  (because  the  functions  H  and  VH  are  singular 
and  highly  oscillatory  as  already  noted),  is  a  trivial  task  within  the  Fourier-Kochin  method.  Indeed, 
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the  Fourier-Kochin  space  integration  consists  in  evaluating  spectrum  functions,  defined  below  by 
(22) ,  given  by  distributions  of  elementary  waves  over  the  surface  EuFuS^.  The  Fourier  integration, 
on  the  other  hand,  is  simpler  in  an  important  respect  but  more  complicated  in  another  respect, 
within  the  Fourier-Kochin  approach  than  the  classical  Green-function  method,  as  is  explained 
further  on. 


The  Fourier  components  in  the  boundary-integral  representations  (14)  are  now  defined  in  terms 
of  Fourier-Kochin  distributions  of  elementary  waves.  The  Fourier  components  0^,  Uy,  (j)^, 
in  the  representation  (14a)  of  time-harmonic  flows  without  forward  speed  can  be  expressed  as 


1 

47r  < 

u^+uy — 
v^A-vy  — 

1  roo  poo 

>  =  lim  —  d/3  da  < 

s  ►-i-O  TT  J — oo  J —oo 

—ia 

y 

-iP 

w^A-  uPf-  — 

<  J  J 

k 

\ 

kDe 


(20a) 


Here,  Dg  =  —k  +  if  e  m  accordance  with  (4) ,  E*  is  given  by 

£*  _  gfeC-i(a^  +  /?i?) 

and  the  function  is  the  wave-spectrum  function  defined  in  the  next  section.  If  F  >  0,  the 
Fourier  components  in  (14b)  and  (14c)  can  be  expressed  as 


<p^+  (j)p  —  (t>^ 


Att  < 


-h  Up  —u^ 

^  > 

+Vp  +vy 

-|-  Wp  -f  4-  ^ 


lim  — 

1-0  TT 


where  Dg  is  given  by  (4) ,  and  S  is  defined  as 


1 

'  sf  ' 

—ia 

qW  ^ 

-ipsi/k 

-ip 

k 

i  a  Si/ k 

k 

0 

\ 

(20b) 


(20c) 


Here,  D  =  {f —Fa)^  —  k  and  the  wave-spectrum  function  is  given  in  the  next  section. 


The  Fourier-Kochin  representation  of  free-surface  effects  given  by  (20)  and  expressions  (21)  and 
(22)  for  the  wave-spectrum  function  is  obtained  in  Noblesse  (2000)  by  substituting  (13b)  into 
(16)-(18) ,  exchanging  the  order  in  which  the  integrations  with  respect  to  the  Fourier  variables 
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(a , /?)  and  the  space  variables  {x  ,y  ,z)  are  performed,  and  performing  several  mathematical  trans¬ 
formations  of  the  resulting  spectrum  functions. 


Wave-spectrum  functions 


The  wave-spectrum  function  5^  in  (20a)  and  (20c)  is  defined  in  Noblesse  (2000)  as 


s  +  st 

O 

II 

'fc _ 

S  +  F^S^ 

>-SP  if< 

f  =  o 

fF>0 

< 

where  the  five  spectrum  functions  S ,  St,  S* ,  S^,  are  given  by 


S  =  j^dA[u  -n  +i^{uxnY-i^{uxnY]e^^  £q 


fs/j 

id){at'^  —  jdt^)/  k 

u-t  ► 

a0/k'^ )  u  ■  t  —  (tyS^u-  V 

f  dApSo 


with  = 


(22a) 

(22b) 

(22c) 

(22d) 


The  spectrum  function  S ,  the  spectrum  functions  St,  S* ,  5^,  and  the  spectrum  function  are 
respectively  given  by  distributions  of  elementary  waves  over  the  boundary  surface  S ,  the  boundary 
curve  r ,  and  the  mean  free  surface  The  spectrum  functions  S ,  are  defined  in  terms  of 

the  normal  components  u-n,  u-V  and  the  tangential  components  u  x  n,  u  -t  of  the  velocity  u  at 
the  boundary  surface  E  and  the  boundary  curve  F.  The  spectrum  function  is  defined  in  terms 
of  the  potential  (f) . 


In  the  special  case  /  =  0 ,  i.e.  for  steady  flow,  the  spectrum  function  is  defined  by  (21)  in 
terms  of  the  spectrum  functions  S  and  S^,  which  do  not  involve  the  potential  (j)  and  are  defined 
directly  in  terms  of  the  velocity  u.  If  />0,  (21)  shows  that  involves  the  spectrum  function 
st  defined  in  terms  of  the  potential  <p ,  although  the  free-surface  boundary  conditions  (3)  can  be 
used  to  express  (f)  in  terms  of  u . 
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Green  functions  and  super  Green  functions 


Expressions  (21)  and  (22)  show  that  the  spectrum  functions  and  S  in  (20)  are  given  by 
distributions  of  elementary  waves  over  the  surface  E  U  F  U  This  surface  can  be  divided  into 
a  set  of  patches  associated  with  reference  points  {Xp  ,yp,Zp  <  0)  located  near  the  centroids  of  the 
patches.  The  spectrum  functions  and  S  in  (20)  are  then  of  the  form 


p=N 

5p  gfc2p+i(Qip+0j/p) 

p=i 

where  summation  is  performed  over  all  the  N  patches  that  represent  the  surface  E  U  F  U  E^,  and 
Sp  is  a  slowly-varying  function  of  a  and  jS .  Expressions  (20a)  and  (20b)  therefore  involve  double 
Fourier  integrals  of  the  type 


e^+Oj-ooJ-oo  Ds 

where  Ap  is  not  a  rapidly  oscillatory  function  of  a  and  /?.  This  Fourier  integral  is  similar  to 
the  Fourier  integral  (11a)  that  defines  the  free-surface  component  in  the  ship-motion  Green 
function.  Indeed,  (11a)  and  (2b)  show  that  is  defined  in  terms  of  the  Fourier  integral 

J—OO 


Expression  (4a)  shows  that  the  foregoing  Fourier  integrals  are  of  the  form 


G  = 


lim 

e-^+O 


(‘oo  yoo  ^  g— z  (  Xa -f  y /3) 

da 


/oo  ro< 

dp  d 

-oo  J—C 


D  +  is 


(23) 


with  A  =  i4p  e  ^  ^  )  or  A  =  A^e  ^  ^  ^ .  A  function  Q  associated  with  a  distribution  of  singu¬ 

larities  is  called  a  super  Green  function  because  of  its  similarities  and  differences  with  usual  Green 
functions,  which  correspond  to  a  point  source. 


The  amplitude  function  A^  corresponding  to  a  point  source  is  simpler  than  the  amplitude  func¬ 
tion  Ap  associated  with  a  Fourier-Kochin  distribution  of  singularities  over  a  patch.  In  fact,  (13c) 
shows  that  the  amplitude  function  corresponding  to  the  Rankine-Fourier  decompositions  (13a)  of 
the  Green  functions  associated  with  the  free-surface  conditions  (3)  is  equal  to  1  if  F  >  0  or  f^k 
if  F  =  0 .  Thus,  the  Fourier  integral  (23)  is,  in  this  respect,  more  complicated  for  a  distribution 
of  singularities  than  for  a  point  source,  i.e.  for  the  Fourier-Kochin  approach  than  for  the  classical 
Green-function  method. 


The  amplitude  function  A  in  (23)  vanishes  exponentially  in  the  limit  fc  oo  if  ^  .Zp  <  0  or 
C  -t-  2  <  0 .  However,  the  amplitude  functions  corresponding  to  a  point  source  and  a  distribution 


16 


of  sources  differ  in  an  important  respect  in  the  limiting  case  C  +  '2p  =  0orC  +  2  =  0.  Specifically, 
the  amplitude  function  A  corresponding  to  a  distribution  of  singularities  over  the  free-surface  plane 
2  =  0  vanishes  as  A:  — >  oo  in  all  cases  (i.e.  even  if  (^  =  0) ,  whereas  ^4  =  1  for  a  point  source 
if  F  >  0  and  C  +  2  =  0 .  This  difference  between  the  Fourier  integrals  associated  with  a  point 
source  and  a  distribution  of  sources  corresponds  to  the  fact  that  the  component  that  accounts  for 
free-surface  effects  is  singular  for  a  point  source  in  the  limit  C  +  2  0  but  is  finite  for  a  distribution 

of  singularities.  Thus,  the  Fourier  integral  (23)  is  in  this  respect  simpler  for  a  distribution  of 
singularities  than  for  a  point  source. 

A  reliable  and  practical  method  for  evaluating  super  Green  functions  of  the  type  defined  by  (23) 
evidently  is  a  critical  requirement  of  the  Fourier-Kochin  representation  of  free-surface  effects.  This 
fundamental  issue  has  been  considered  for  a  generic  dispersion  function  D  +  i  sD'  and  a  generic 
amplitude  function  A,  i.e.  for  a  broad  class  of  dispersive  waves  and  arbitrary  distributions  of  sin¬ 
gularities,  in  several  studies.  The  most  important  result  of  these  studies  of  generic  super  Green 
functions  —  obtained  by  the  author  in  cooperation  with  Dr.  Chi  Yang  and  first  given  in  Noblesse 
and  Chen  (1995)  —  is  a  formal  decomposition  of  the  generic  super  Green  function  (23)  into  wave 
and  local  components,  and  a  remarkably  simple  Fourier  representation  of  the  wave  component. 
Alternative  Fourier  representations  of  the  local  component  in  this  waved-local  decomposition  are 
given  in  Noblesse  and  Yang  (1996) ,  Noblesse  and  Chen  (1997) ,  and  Noblesse  et  al.  (1999) .  The 
most  practical  of  these  alternative  representations  is  that  given  in  Noblesse  and  Chen  (1997) .  This 
study  and  Noblesse  et  al.  (1999)  are  summarized  below. 

Practical  representation  of  generic  super  Green  functions 

The  double  Fourier  integral  (23)  has  been  considered  in  numerous  studies  of  free-surface  Green 
functions.  These  studies  rely  on  integration  contours  in  the  complex  Fourier  plane  a  =  ar  +  iai,  or 
alternatively  in  the  complex  Fourier  planes  /3  =  Pr  +  i or  k  =  kr  +  ih ,  in  the  manner  explained 
in  Noblesse  (1981) .  A  different  approach  is  used  here  and  in  the  previously-listed  studies  of  generic 
super  Green  functions  :  the  Fourier  representation  (23)  is  analyzed  directly,  i.e.  without  introducing 
a  contour  of  integration,  in  the  real  Fourier  space  {a,P) . 

The  double  Fourier  integral  (23)  can  be  expressed  in  the  basic  form 


Q  =  G2  +  Gi 


(24a) 
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Here,  Q2  is  the  double  Fourier  integral 

=  (24b) 

J—ooJ—oo 

obtained  by  formally  setting  e  =  0  in  (23) .  Noblesse  et  al.  (1999)  show  that  Qi  =  Q  -  Q2  can  be 
expressed  in  terms  of  single  Fourier  integrals  along  the  dispersion  curve(s)  defined  by  the  dispersion 
relation  D{a ,  /3)  =  0 : 

=  «  (240 

II  VD||  =  yjDl+Dl  where  and  Dj}  are  the  derivatives  of  the  dispersion  function  D  with  respect 
to  the  Fourier  variables  a  and  fi ,  and  ds  stands  for  the  differential  element  of  arc  length  along  a 
dispersion  curve. 

A  fundamental  result  of  the  theory  of  Fourier  transforms,  given  in  Lighthill  (1970) ,  is  that  the 
behavior  of  a  function  f{x)  in  the  farfield  a:  — >  00  is  determined  by  the  dominant  singularities  of 
its  Fourier  transform.  Thus,  the  dominant  farfield  features  of  the  function  Q2  defined  by  (24b) , 
i.e.  the  farfield  waves  contained  in  Q2 ,  stem  from  the  dispersion  curves  defined  by  the  dispersion 
relation  D  =  0 .  Specifically,  Noblesse  et  al.  (1999)  yield 


Q2 


f ds^^A 
ko  ■'o-o 


,-i{Xa  +  Yp) 


(25a) 


in  the  farfield  limit  +  Y"^  00 .  Here,  erf  is  the  usual  error  function  and  ip  is  defined  below. 

The  basic  decomposition  (24a)  and  the  farfield  approximation  (25a)  yield 


G-^G^=Gi+GV 


(25b) 


as  x2+y2 

— »  00 .  Expressions  (24c)  and  (25a)  yield 


sign(£>0  +  erf(^)  i^xa  +  Yp) 
||VZ)|| 


The  function  ip  in  (26a)  is  given  by 


,  ^  k{XDa+YD/3) 

^  CTllVPIl 

where  cr  w  1  is  a  positive  real  number.  The  error  function  erf(V’)  becomes 


(26a) 


(26b) 


eTi{ip)r^  sign{X Da+YD/s)  as  X^+Y^^oo 


(26c) 
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The  term  sign(T>')  in  (26a) ,  which  stems  from  the  component  Qi  in  (24a) ,  ensures  that  the 
radiation  condition  is  satisfied. 


In  the  nearfield,  the  super  Green  function  Q  can  be  expressed  as  the  sum  of  the  wave  component 
given  by  (26)  and  a  local  component  : 

G  =  g^+g^  (27) 

Expressions  (24a)  and  (25b)  show  that  the  local  component  G^  =  G  —  G^  is  given  by 


Alternative  representations  of  the  local  component  g^  in  the  fundamental  decomposition  (27)  into 
wave  and  local  components  are  given  in  Noblesse  and  Yang  (1996) ,  Noblesse  and  Chen  (1997) , 
and  Noblesse  et  al.  (1999) .  The  latter  study  considers  special  cases,  corresponding  to  simple 
open  and  closed  dispersion  curves,  and  the  general  case  of  arbitrary  dispersion  curves.  A  simpler 
representation  that  is  valid  for  arbitrary  dispersion  curves  is  (see  Appendix  3) 

g-i(Xa  +  y/3) 

D=0 


TOO  roo 

g’-=  dn  ho  (A -Y,  El  Aj) 

J—OOJ—OO  r,  „ 


D 


(28a) 


Summation  in  (28a)  is  performed  over  all  the  dispersion  curves,  which  are  identified  by  the  subscript 
j ,  and  Ej  is  defined  as 


\\VD\\^^^^Uk]\\VD\\] 


DDi 


+  iA(a-a,-|^)  +  *F(/3-/3,.- 


)] 


(28b) 


The  subscript  or  superscript  j  in  (28)  indicates  evaluation  at  the  dispersion  curve;  e.g.  Aj 
stands  for  the  value  of  the  amplitude  function  A  at  the  dispersion  curve. 


The  local  component  g^  defined  by  (28)  can  be  further  decomposed  as  in  Noblesse  and  Chen 
(1997)  : 


G^=G^+gh 

(29a) 

Here,  the  component  g^  is  defined  as 

roo  poo  1  Z? 

gi:=  fdp  fda^  ^  + 

J—ooJ—oo  D 

(29b) 

w.th  «  = 

(29c) 

The  component  g^  in  (29a)  is  then  given  by 

gh  =  [dl3  g-i(Xa  +  Y0) 

J  ~oo  J—OO 

(29d) 
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where  the  functions  Ej  and  E  are  given  by  (28b)  and  (29c) .  The  component  defined  by  (29b) 
accounts  for  the  contribution  of  the  entire  Fourier  plane  except  the  regions  near  the  dispersion 
curves.  The  contribution  of  these  dispersion  strips  is  given  by  the  component  defined  by  (29d) . 
The  integrand  of  the  double  Fourier  integral  (29b)  vanishes  at  a  dispersion  curve  D  =  0,  and  this 
Fourier  integral  can  be  evaluated  simply  and  very  efficiently  if  the  amplitude  function  A  is  not 
rapidly  oscillatory,  as  is  assumed  here.  The  integrand  of  (29d)  is  finite  at  a  dispersion  curve,  and 
the  component  Qp  can  also  be  evaluated  without  difficulty. 


In  fact.  Noblesse  and  Chen  (1997)  shows  that  this  component  can  be  approximated  by  the  single 
Fourier  integrals 


2y/n 


—  Vfds 


e-^(^a  +  y0) 

livz)rai/  > 


(30a) 


along  the  dispersion  curves  D  —  0.  Here,  ^  is  given  by  (26b)  and 


C  =  {DaaDl  +  2Da0D^D0  +  )/||VL>f  (30b) 


Furthermore,  dA/dv  is  the  derivative  of  the  amplitude  function  A{a ,  0)  in  the  direction  of  the  unit 
vector  u  =  VD/||VD1|  normal  to  a  dispersion  curve  in  the  Fourier  plane  {a,0) .  The  derivative 
dA/dv  can  be  determined  using  analytical  or  numerical  differentiation,  i.e. 

dA  ^  DaAaAD^Ap  _  A{a\0+)-A{a-,0-) 
dv  ||VD||  2  b 

Here,  A(a*,/3^)  are  the  values  of  the  function  A  at  points 

(a=t,^±)  =  (a,^)±6(Z)a,r>/3)/||VD|| 

defined  via  displacements  ±bv  from  a  point  [a  ,0)  of  a  dispersion  curve. 

In  summary,  (27)  and  (29a)  show  that  the  super  Green  function  Q  defined  by  the  singular  double 
Fourier  integral  (23)  can  be  expressed  as 

g  =  g^+ =  (3i) 


(30c) 


(30d) 


g^  and  g^  are  given  by  the  single  Fourier  integrals  (26)  and  (30) ,  and  is  defined  by  the  regular 
double  Fourier  integral  (29b)  and  (29c) .  The  single  Fourier  integrals  (26a)  and  (30a)  can  evidently 
be  combined,  as  in  Noblesse  and  Chen  (1997) .  The  wave  component  ^^and  the  local  components 
g^  and  g^  in  the  representation  (31)  can  be  evaluated  in  a  straightforward  and  highly  efficient 
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manner  if  the  amplitude  function  A  in  (23)  is  not  rapidly  oscillatory. 


Fourier-Kochin  representation  of  nearfield  waves 

The  decomposition  (27)  shows  that  the  double  Fourier  integrals  (20a)  and  (20b) ,  which  are  of 
the  form  (23) ,  can  be  expressed  as  sums  of  wave  and  local  components  ; 

(32a) 


(32b) 


-  (^P  1 

u^+  ay  —  np  j 

[u^+u^ 

] 

^  ^ 

nPj 

[a^+u\ 

(l>^+<f>F  +(l>r  -4'^ 
u^+u^  +  u^+up 

The  wave  components  in  the  decompositions  (32)  are  defined  by  (26a)  in  terms  of  single  Fourier 
integrals  along  the  dispersion  curves  defined  by  the  dispersion  relation  D  =  0 .  It  follows  that  the 
second  component  on  the  right  side  of  (20c)  does  not  contribute  to  the  waves  contained  in  (20b) . 

Specifically,  the  wave  component  (j)^ ,  in  (32)  is  given  by 


_ 

( 

1 

47r  < 

w 

w 

V 

>=-<e/ 

^qJd=o 

ds 

IlVDII  ' 

—iot 

-ip 

w 
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El 
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c* 

■'p 


(33a) 


The  identity  /^  =  fc  at  a  dispersion  curve  D  =  0 ,  which  follows  from  (5a) ,  was  used  in  (20a) . 
The  functions  D  and  D'  in  (33a)  are  given  by  (4b) .  Summation  is  performed  over  all  the  patches 
that  represent  the  surface  S  U  F  U  as  previously  explained.  The  function  stands  for  the 
contribution  of  patch  p  to  the  wave-spectrum  function  defined  by  (21)  and  (22) ,  with  the 
substitutions 

(s,r,E^)-^(Sp,rp,s;) 

(^x ,y , z)  >  (x  Xp ,y  yp Zp) 
performed  in  (22) .  The  function  £*  is  given  by 


(33b) 


5*  =  e(C  +  zp)fc-i[(4-®p)a-i-(»?-yp)/3] 


(33c) 


Finally,  ipp  is  defined  by  (26b)  as 


/  —  ]c  ( C  )  -Pq  d-  ( ?7  Vp)  DjS 


(33d) 
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In  the  farfield,  we  have 


erf(V'p)  sign[ -  Xp)  Da+ {r}  -  yp)  Dp] 

and  division  of  the  surface  E  U  F  U  into  patches  may  not  be  necessary,  as  shown  in  Yang  et 
al.  (2000a) .  In  summary,  farfield  and  nearfield  waves  are  defined  by  (33)  in  terms  of  the  wave- 
spectrum  function  5^^  given  by  (21)  and  (22) . 


Time-harmonic  ship  waves 


The  dispersion  curves  in  the  regions  (7)  yield  distinct  contributions  to  the  wave  component  , 
defined  by  (33a) .  If  r  <  1/4 ,  (j)^  and  can  be  expressed  as 

M 

The  components  respectively  correspond  to  the  dispersion  curves  in  the  inner  region 

Oi~  <Q<ocf  and  the  outer  regions  a<a~  and  Q+<a .  The  component  represents  a  system  of 
ring-like  waves,  called  ring  waves.  The  components  and  ujy  represent  two  systems  of  V-like 
waves  contained  within  wedges  and  called  outer-V  waves  and  inner-V  waves,  respectively.  Curves 
of  constant  phase  (e.g.  crest  lines)  corresponding  to  the  ring  waves  (thick  solid  lines),  the  outer-V 
waves  (dashed  lines),  and  the  inner-V  waves  (thin  solid  lines)  are  depicted  on  the  left  side  of  Fig. 
3  for  t  =  0.24 .  In  the  limit  t  =  0,  the  components  and  respectively  correspond  to 

time-harmonic  ring  waves  (limit  F=0)  and  steady  ship  waves  (limit  /  =  0) . 


<  +  Cv+<?K' 


w 


kW 


u 


+  u^+ul 


^oV 


W 


iV 


(34a) 


If  r  >  1/4 ,  (f>^  and  can  be  expressed  as 


_ y 

>  =  i 

(34b) 


where  u^p  and  are  associated  with  the  dispersion  curves  in  the  regions  a<a'l  and  a'^  <a, 
respectively.  The  component  represents  a  system  of  inner-V  waves  qualitatively  similar  to  the 
inner-V  waves  in  (34a) .  The  component  u^p  represents  fan-like  waves  and  incomplete  ring  waves. 
These  waves  are  also  contained  within  a  wedge,  and  are  called  ring-fan  waves.  The  system  of 
ring-fan  waves  can  be  further  divided  into  two  systems  of  waves  :  a  system  of  inner-fan  waves  and 
a  system  of  partial-ring  waves  and  outer-fan  waves,  which  correspond  to  the  portions  o;  <  —f/F 
and  —f/F<a<af  of  the  dispersion  curve  a<af.  Constant-phase  curves  corresponding  to  the 
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partial-ring  waves  and  the  outer-fan  waves  (thick  solid  lines),  the  inner- fan  waves  (dashed  lines), 
and  the  inner-V  waves  (thin  solid  lines)  are  depicted  on  the  right  side  of  Fig.  3  for  t  =  0.26  . 

Expressions  for  the  four  basic  time-harmonic  ship-wave  components  —  i.e.  the  inner-V  waves  in 
(34a)  and  (34b) ,  the  ring-fan  waves  in  (34b) ,  and  the  outer-V  waves  and  the  ring  waves  in  (34a)  — 
can  be  obtained  from  (33)  and  are  now  given.  The  inequalities  (8)  yield  sign(I?')  =  sign(/— Fa)  =  1 
for  all  the  dispersion  curves  located  in  the  regions  defined  by  (7) ,  except  the  dispersion  curve  in 
the  region  a+  <  a  for  which  sign(F>')  =  —1 .  The  ship-motion  dispersion  function,  defined  by  (4b) 
as  D  =  {Fa  —f)^  —  k ,  yields  =  —0/k  and  Da  =  —a/k  with 

a  =  a  -t-  2  fc  (r  —F^a) 


Expression  (33d)  then  yields  ipp  =  —ipp  with 


..  _  7-  i^-Xp)a  +  {r]-yp)0 

Pp  -  K  /■-" . 


(35) 


The  dispersion  curves  in  the  three  regions  a+  <  a,  a  <  af ,  a  <  a~  can  be  defined  by  the 
parametric  representation 

[  ^  ] 

a  =  fjF-eVklF 


/3  =  Vk^  —a^  sign(f) 


where  -oo<t<oo  and 
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with  a+  ,  a^  ,  a^  given  by  (6) .  a  in  (35)  is  then  given  by 

a  =  f/F  +  e  {2F^k -l)Vk/F 

Finally,  (33a)  and  the  relation  ds/||VDl|  =  sda/Dp  =  —ekda/0  yield 
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(36a) 


(36b) 


(36c) 
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(36d) 


The  subscript  7  in  and  stands  for  iV  for  inner-V  waves,  RF  for  ring-fan  waves,  or  oV  for 


outer  V-waves.  (36a)  yields  (3  ~ 

tjsfC  as  t— >0  with 

1+1/\/T+47 

inner-V 

C  =  < 

1  — l/'\/l+  4r 

^  for  < 

ring-fan 

l  +  l/'s/l—  4r 

outer-V 

s.  > 

>  waves 


Thus,  the  integrand  of  (36d)  is  continuous  at  t  =  0. 

The  dispersion  curve  in  the  region  aj  <a<af  can  be  defined  by  the  parametric  representation 


(37a) 


k  =  /^/  (v'  1/4  ■+  Tcost  +1/2)2 
a  =  k  cost  (3  =  k  sint 
with  -7r  <  f  <  TT .  In  the  limit  r  ^  0 ,  the  inner  dispersion  curve  becomes  a  circle  of  radius  k  =  f^ 
centered  at  the  origin  of  the  Fourier  plane.  3  in  (35)  is  given  by 

a  =  k[2T  +  {l-2F^k)  cost] 

(33a)  and  the  relation  ds/||VD||  =  —kdt/Dk  yield 


(37b) 
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(37c) 


The  integrand  of  (37c)  is  continuous  for  “7r<i <7r  since  ring  waves  exist  only  for  r  <  1/4 .  In  the 
limit  F  =  0,  expressions  (37)  for  the  ring  waves  become 
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(38) 


with  {k,a,P)  =  /2(1,  cost ,  sint )  and  (/jp  =  [  (^  -  Xp )  cost  +  {r]  -  yp)  smt]f‘^/a .  (38)  defines 
time-harmonic  waves  associated  with  diffraction-radiation  without  forward  speed. 

Both  the  inner-V  waves  and  the  outer-V  waves  in  (34a)  contain  real  and  imaginary  parts. 
However,  the  waves  given  by  the  sum 
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are  real  in  the  special  case  /  =  0 ,  i.e.  for  steady  ship  waves.  Specifically,  (36)  and  the  property  that 
steady-flow  spectrum  functions  satisfy  the  relation  S^{-a ,  -/?)  =  S^{a ,  /?) ,  with  Z  =  complex 
conjugate  of  Z ,  show  that  steady  ship  waves  can  be  expressed  as 

User^ 

yW  J-oO  k-V 

w 

J 

3Je  stands  for  the  real  part,  the  relation  ds/||VD||  =  d^/Da  =  —kdfd/a  was  used,  and 

k  =  v+  \jv^+  a  =  v/fc/F  (39b) 

with  V  =1/(2F^) .  a  in  (35)  is  given  by  3  =  (1—  2F^k)a. 

Rankine-Fourier  decomposition  :  physical-space  analysis 
Consider  the  free-surface  boundary  condition  (9) .  Define  the  coordinates  (X,  Y,  Z)  as 

{i-x,r)-y,(:  +  z)  =  K{X,Y,Z) 

where  A  stands  for  a  length  scale.  With  respect  to  the  A-scaled  coordinates  (X,  Y,  Z) ,  the  free- 
surface  condition  (9)  becomes 

Gz-fKG  +  i2rGx  +  {F‘^lk)Gxx  =  Q  (40) 

In  the  particular  case  F=0 ,  (40)  becomes  Gz  —  pKG  =  0 .  This  condition  yields 

Gz^o]  f  A  <  1//2 1  f  /  <  1//2 

^  if  <  >  i.e.  for 

G  «  0  A  »  1//2  /  >  1//2 

Thus,  for  time-harmonic  flow  without  forward  speed,  the  free  surface  acts  like  an  infinitely  rigid 
wall  in  the  nearfield  region  r'  l//^  and  like  an  infinitely  soft  wall  in  the  farfield  region  /  »  1//2. 
The  alternative  representations  (12)  therefore  are  complementary  and  appropriate  in  the  nearfield 
and  the  farfield  : 

r'  <  1//2 
r'  >  1//2 

These  nearfield  and  farfield  approximations  indicate  that  an  ’’optimal”  Rankine-Fourier  decompo¬ 
sition  of  the  Green  function  G  related  to  the  free-surface  boundary  condition  (3c)  is 

47r  G  =  -1/r  -  1/r'  -H  2/r"  +  T~  (41a) 


-l/r-l/r'  +  7f 
-1/r -h  1/r' +  ff 


25 


where  r'  and  r"  are  defined  by  (lb)  with  (7^=  1//^,  and  T  is  a  Fourier  component. 

In  the  particular  case  f  —  0,  (40)  becomes  Gxx  +  (A/F^)  Gz  —  0 .  This  condition  yields 


Gxx  ^  0  ] 

>  if  ^ 

f  A  <  1 

i.e.  for  ^ 

f  r'  <  F2 

Gz^O  j 

A  »  F2  J 

r’:$>F^ 

<  j 

Thus,  for  steady  flow,  the  free  surface  acts  like  an  infinitely  soft  wall  in  the  nearfield  region  r' 
and  like  an  infinitely  rigid  wall  in  the  farfield  region  r'  >  F'^.  The  alternative  representations  (12) 
therefore  are  appropriate  in  the  nearfield  and  the  farfield  : 


47rG=  I 


-1/r  -H  1/r'  -h  H  ] 

y  for 

1  r'  <  F2  1 

-l/r-l/r'  +  n 

<  > 

|r'»F2j 

These  nearfield  and  farfield  approximations  indicate  that  an  ’’optimal”  Rankine-Fourier  decompo¬ 
sition  of  the  Green  function  G  related  to  the  free-surface  condition  (3b)  is 

47r  G  = —1/r -j- 1/r' —  2/r"  +  (41b) 

where  r'  and  r"  are  defined  by  (lb)  with  C‘^  =  F'^,  and  is  a  Fourier  component. 

If  fF  >  0 ,  the  free-surface  condition  (40)  yields 


Gxx^O 

[a^oI 

'  r'-^O  1 

^  if  < 

1  j 

>  i.e.  for  < 

1  / 

GsaO  J 

[A^ooJ 

I  r'— >  00  1 

Thus,  the  free  surface  acts  like  an  infinitely  soft  wall  in  both  the  nearfield  r'  — +  0  and  the  farfield 


r'  — >  00 ,  and  the  Rankine-Fourier  decomposition  of  the  Green  function 

47rG  =  — 1/r  +  l/r'  +  H 

is  ’’optimal”  for  the  ship-motion  free-surface  condition  (3a)  with  fF  >  0. 
Expressions  (41)  yield 


(41c) 


’-l/r'+2/r"' 

f  ■> 

F  =  0 

47rG  =  —1/r  -|-^ 

l/r'-2/r" 

>  +  F 

if  < 

/  =  0  ^ 

1/r' 

<  > 

/F>0 

V. 

(42a) 


Here,  r ,  r',  r"  are  defined  by  (lb)  with  G^  given  by 


C2  =  < 

[f^o 
mM  ? 

V=o| 

(42b) 
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The  Fourier-component  T  in  (42a)  is  given  in  the  next  section. 


Rankine-Fourier  decomposition  :  Fourier-space  analysis 

The  free-surface  component  in  expression  (10a)  for  the  Green  function  G  corresponding  to 
the  free-surface  boundary  condition  (9)  is  given  by  (11) .  Alternative  forms  of  (11)  are 


47rG^  =  1  +  F+  =  4-HF-  =  A-  4+:^+  =  -4  +  4+:F- 


(43) 


y**  ytii 

The  functions  G^,  T'^  are  given  by  the  right  side  of  (13b)  with  k,  £ ,  Dg  given  by  (2b)  and 

(4) ,  and 


A  =  l 


l(f-Faf/k+l]/2 

A+ 

1 

A- 

^  < 

U-Faf/k  > 

A+ 

l+[{f-Fotf/k-l]e-<^^'^ 

A- 
<  > 

{f-Fa)‘^/k-[{f-Faf/k-l]e-<^^^ 

(44a) 


as  follows  from  (11),  (2c)  and  (2d) .  The  functions  H+  and  H~  in  (43)  are  identical  to  the 
previously-defined  functions  H  and  Ti.  (the  notation  is  used  here  for  convenience) . 


The  theory  of  Fourier  transforms,  Lighthill  (1970) ,  shows  that  the  farfield  behaviors  of  the 
local  components  contained  in  the  functions  G^,  are  determined  by  the  behaviors  of  the 

functions  A^/D ,  A^/ D ,  D  at  the  origin  A:  =  0  of  the  Fourier  plane.  Conversely,  the  behaviors 
of  these  functions  in  the  limit  k  oo  determine  the  nearfield  behaviors  of  the  functions  G^,  H^, 
in  the  physical  space. 


In  the  particular  case  F  =  0,  we  have  C^  =  l//2  and  (44a)  becomes 


f  > 

A^ 

(/VA;+l)/2 

A+ 

1 

A- 

^  =  < 

f/k 

A+ 

1+  (f^/  k  —1)  e~^  ! 

A- 

> 

k  -  (/^/A;  -l)e“''/'^^ 

(44b) 


The  functions  A^,  A'^  and  A'^  are  finite  in  the  limit  A:  — >  oo ,  whereas  the  functions  A~  and  A~  are 
C*(/^/  k) ,  i.e.  vanish,  in  this  limit.  It  follows  that  the  functions  H~  and  T~  are  less  singular  at  the 
origin  r'  =  0  than  the  functions  G^,  and  JF+  in  (43) .  In  the  limit  A:  0 ,  the  function  A~  is 
0{Plk) ,  i.e.  singular,  whereas  A~  is  0(1) ,  i.e.  finite.  It  follows  that  the  local  component  contained 
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in  the  function  T  decays  faster  in  the  farfield  oo  than  the  local  component  associated  with 
the  function  }i~  .  The  decomposition  (41a)  is  therefore  preferable. 

In  the  particular  case  /  =  0 ,  we  have  (7^  =  and  (44a)  becomes 


(FV/A:+1)/2 

4+ 

1 

A- 

^  ^ 

F^a'^/k 

A+ 

l+(FV/A:-l)e-'^"'' 

A- 

<  > 

F V/ -  ( F  V/ A;  - 1 )  e-'P’"  ^ 

k.  / 

In  the  limit  oo ,  the  functions  A^,  A~  and  A~  are  0{F‘^a^/  k) ,  i.e.  unbounded,  whereas  the 
functions  A'^  and  are  0(1) .  The  functions  and  therefore  are  less  singular  at  the  origin 
r'  =  0  than  the  functions  H~  and  in  (43) .  In  the  limit  k  ^  0,  the  function  A'^  is  finite 
but  A'^  is  0{F^k)  and  thus  vanishes.  The  local  component  contained  in  the  function  F"*"  therefore 
decays  faster  in  the  farfield  r'  — >  oo  than  the  local  component  corresponding  to  the  function  H'^  , 
and  the  decomposition  (41b)  is  preferable. 

In  the  general  case  /F  >  0 ,  the  functions  A^,  A~  and  A~  in  (44a)  are  0{F‘^a‘^/  k) ,  i.e. 
unbounded,  and  the  functions  A"^  and  A"^  are  0(1)  in  the  limit  k  oo.  In  the  limit  A:  — >  0,  the 
function  >1+  is  finite  but  the  function  .4"^  is  0{f‘^/k),  i.e.  unbounded.  Thus,  the  decomposition 
(41c)  is  preferable  in  this  case. 


The  component  F  in  expressions  (42a)  for  the  Green  functions  related  to  the  free-surface  con 
ditions  (3)  is  then  defined  by  the  Fourier  superposition  of  elementary  waves  (13b) ,  i.e. 

1  f°° 


F=  lim 

e— ^H-0  TT  j—oo  J—oo 


(45) 


Here,  k  and  £  are  given  by  (2b) ,  the  dispersion  function  Dg  is  given  by  (4) ,  and  the  amplitude 
function  A  is  defined  by  (44)  as 


(46) 


'fyk-{fyk-i)e-'^/f^' 

r 

F  =  0 

A=  < 

l+(FV/A:-l)e-'^''' 

>  if  < 

/  =  0  > 

1 

k.  j 

/F>0 

Expressions  (46)  yield 


< 


(H 

>  as  k/f"^^  < 

foj 

l/VtJ 

[ocj 

S  if  F  =  0 
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{F‘^k{l+a^lk‘^)\ 

I  1 


>  as  F^k-^{ 


0 

00  I 


if  /  =  0 


In  the  general  case  fF>0  and  the  special  cases  /=0  and  F  =  (46)  yield 


A 

D 


'  0(1)  I 

'o(l/fc2)J 


as  k 


and  >1  =  1  if  Z)  =  0.  Thus,  the  amplitude  functions  A  defined  by  (46)  are  equal  to  1  at  a  dispersion 
curve  D  =  0,  and  the  free-surface  components  defined  by  (46)  and  (13)  yield  identical  farfield  waves 
as  one  expects. 


Expressions  (42) ,  (45)  and  (46)  define  the  Green  functions  associated  with  the  free-surface 
boundary  conditions  (3)  in  terms  of  elementary  Rankine  singularities  and  Fourier  superpositions 
of  elementary  waves.  The  Rankine-singularity  components  in  (42)  account  for  the  dominant  terms 
in  both  the  nearfield  and  farfield  asymptotic  approximations  to  the  nonoscillatory  local  compo¬ 
nents  included  in  the  Green  functions  corresponding  to  (3) .  The  foregoing  method  of  extracting 
the  leading  nearfield  and  farfield  local  component  from  the  Fourier  component  can  in  principle  be 
pursued  to  obtain  additional  terms  in  the  nearfield  and  farfield  asymptotic  approximations  to  the 
local  components  in  free-surface  Green  functions. 


Rankine  and  Fourier-Kochin  representation  of  nearfield  flows 

The  Rankine-Fourier  decomposition  process  that  has  been  expounded  above  for  the  simple  case 
of  free-surface  Green  functions  is  useful  more  widely,  notably  within  the  Fourier-Kochin  approach. 
In  fact,  the  potential  and  velocity  representations,  the  Fourier-Kochin  approach,  and  the  Fourier- 
Rankine  decomposition  process  yield  remarkably  simple  analytical  representations  of  nearfield  flows 
for  the  three  basic  types  of  free-surface  fiows  corresponding  to  (3) ,  as  is  now  explained. 

Transformations  of  the  wave-spectrum  functions  5^,  guided  by  the  Rankine-Fourier  decom¬ 
position  process,  can  be  performed  for  the  purpose  of  obtaining  practical  representations  of  the 
local  component  in  (32) .  These  transformations  express  the  wave-spectrum  functions 

given  by  (21)  and  (22)  in  terms  of  local-spectrum  functions  and  components  that  involve  the 
dispersion  function  D  and  are  dominant  in  both  the  limit  k  ^  oo  and  the  limit  fc  ^  0 : 

S^=DS^+S^  (47) 

The  dominant  component  D  can  be  expressed  in  terms  of  distributions  of  Rankine  singularities. 
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The  potential  (j)  and  velocity  u  defined  by  the  boundary-integral  representations  (14)  and  the 
Fourier-Kochin  representations  (32)  of  free-surface  effects  can  then  be  expressed  as 


/^1 

>  =  ^ 

>  +  ^ 

\u\ 

i“  J 

kR 


/kR 


with 


= 


V 

J 

V 

'J 

(>1 

>  if  ^ 

1  F=0 

l-‘J 

F>0 

k  J 

>  +  < 

1 

P^Up  -\-Uy^ 

f  ^ 

1 

o 

II 

_ 

II 

-1 

►  if  < 

/=0  ^ 

0 

/F>0 

(48a) 


(48b) 


The  component  on  the  right  side  of  (48a)  corresponds  to  the  nearfield  waves  considered 

previously,  and  the  other  components  represent  a  nonoscillatory  local  flow  disturbance.  Thus,  the 
wave  component  u'^is  given  by  single  Fourier  integrals  along  the  dispersion  curves  D  =  0 .  The 
Fourier  integrals  (33a)  involve  the  wave-spectrum  functions  Sidelined  by  (21)  and  (22)  in  terms 
of  distributions  of  elementary  waves  over  the  surface  S  U  F  U 

The  four  components  in  the  second  group  on  the  right  side  of  (48a)  are  defined  in  terms  of 
distributions  of  Rankine  singularities.  Specifically,  the  Rankine  component  </>^,  is  given  by 
(15a)  in  terms  of  distributions  of  elementary  Rankine  singularities  R  =  —1/r  over  the  boundary 
surface  S.  The  Rankine  component  4>^ ,  is  also  given  by  distributions  of  elementary  Rankine 
singularities  over  the  boundary  surface  S;  specifically,  the  component  (f)^ ,  is  given  by  (15b) 
with  R*  defined  as 


R*  =  < 


-l/r'-b2/r" 

< 

^  if  . 

o 

II 

[-1//  J 

[fF>0j 

(49) 


Here,  r'  and  r"  are  given  by  (lb)  with  defined  by  (42b) .  The  function  R*  is  evidently  related  to 
the  Rankine-Fourier  decompositions  (42)  of  the  Green  functions  corresponding  to  the  free-surface 
boundary  conditions  (3) .  The  Rankine  component  (f)^ ,  Up  in  (48a)  stems  from  the  forcing  term  p 
in  (3)  and  is  given  by 


47r  < 


'  ct>f 

'  (1/0" ' 

Up 

< 

►  =  2  /  dA  p  i 

(1/0^' 

Vp 

(i/o;' 

Wp 

J 

.(1/0^ 

The  Rankine  component  p ,  Up  in  (48a)  is  given  by  distributions  of  Rankine  singularities  over 
the  boundary  curve  F . 
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Finally,  the  two  components  in  the  third  group  on  the  right  side  of  (48a)  are  Fourier  components 
that  are  defined  by  single  Fourier  integrals  along  the  dispersion  curves  and  double  Fourier  integrals, 
as  given  by  (30)  and  (29b)-(29c) .  These  Fourier  integrals  are  defined  in  terms  of  the  local-spectrum 
function  in  (47) .  The  integrands  of  the  double  Fourier  integrals  defining  the  component  <j)^  , 
are  finite  at  fc  =  0  and  continuous  ever3rwhere  in  the  Fourier  plane,  and  decay  rapidly  as  fc  — »  oo . 

The  Rankine  component  (/>p,  Up  and  the  Fourier  components  and  (j)^ ,  in  the 

nearfield  flow  representation  (48) ,  called  Rankine  and  Fourier-Kochin  representation,  will  be  given 
elsewhere.  The  Rankine  and  Fourier-Kochin  representation  of  nearfield  flows  does  not  involve  the 
special  functions  that  appear  in  the  integral-representations  of  the  Green  functions  associated  with 
the  free-surface  boundary  conditions  for  steady  and  time-harmonic  ship  waves  given  in  Noblesse 
(1981)  and  Chen  (1999).  In  fact,  the  Rankine  and  Fourier-Kochin  representation  (48a)  is  entirely 
defined  in  terms  of  the  elementary  solutions  (1) ,  and  thus  is  remarkably  simple. 

Iterative  solution  procedures  and  related  approximations 

Numerical  calculations  of  free-surface  potential  flows  due  to  ships  and  offshore  structures  have 
traditionally  relied  on  the  computation  of  influence  coefficients  associated  with  a  discretization  of 
an  integral  equation  obtained  from  the  classical  potential  representation.  An  alternative  approach 
based  on  an  iterative  solution  procedure  can  be  advantageous,  and  can  provide  useful  analytical 
approximations.  An  example  is  the  slender-ship  approximation  given  in  Noblesse  (1983)  and  No¬ 
blesse  and  Triantafyllou  (1983) .  This  approximation  defines  the  steady  velocity  field  u{^)  due  to  a 
ship  explicitly  in  terms  of  the  speed  and  shape  of  the  ship  : 

ff(|)  =  l^dAix)  rf  (x)  V^G({,  x)  -H  F‘^j^dC{x)  [  rf  (f )  f  ty{x)  x)  (50) 

The  slender-ship  approximation  (50)  is  a  generalization  of  the  Michell  thin-ship  approximation. 
Specifically,  the  Michell  approximation  differs  from  (50)  in  that  it  defines  u  in  terms  of  a  distribution 
of  sources,  with  strength  2  n®,  over  the  ship  centerplane  y  =  0  instead  of  a  distribution  of  sources 
of  strength  n®  over  the  port  and  starboard  sides  of  the  actual  ship  surface  E .  In  addition,  there  is 
no  distribution  of  sources  around  the  waterline  F  in  Michell’s  approximation. 

The  wave  component  in  the  decomposition  u  =  is  given  by  the  Fourier-Kochin 

representation  (33) ,  which  can  be  expressed  as  in  Yang  et  al.  (2000b) .  The  wave-spectrum  function 
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5^  corresponding  to  the  slender-ship  approximation  (50)  is  defined  by 

S^  =  J^dAn^e'^^  So  + 

with  So  =  gi(Q*  +  /?2/)  The  local  component  can  be  effectively  evaluated  from  (50)  with  the 
Green  function  G  taken  as  the  simple  analytical  approximation  to  the  local  component  defined 
by  (42)  as 

47rG^  ss  — 1/r  +  1//  —  2/r" 

Here,  r,  r',  r"  are  defined  by  (lb)  with  C^  —  F^. 

At  the  ship  hull  surface  E ,  the  velocity  field  u  given  by  the  slender-ship  approximation  (50)  can 
be  modified  using  the  transformation 


Ut  =  u  —  {u  •  n  —  n^)n  (51) 

This  transformation  yields  u*  -  n  =  and  thus  ensures  that  the  velocity  distribution  u*  at  the 
hull  surface  satisfies  the  hull  boundary  condition.  The  hull-condition  transformation  (51)  can  be 
applied  to  any  velocity  distribution  u  computed  at  a  ship  hull  surface.  In  particular,  (51)  shows 
that  the  velocity  distribution  u*  associated  with  the  trivial  velocity  distribution  u  =  0  is  given  by 
u*  =  n^n.  This  velocity  distribution  is  normal  to  the  ship  surface  E  and  evidently  satisfies  the 
hull  boundary  condition  u*  -  n  =  n®.  The  velocity  distribution  «»  =  n®  n  in  fact  corresponds  to  the 
slender-ship  approximation  (50) .  The  practical  usefulness,  notably  for  hull-form  optimization,  of 
this  remarkably  simple  approximation  is  demonstrated  in  Letcher  et  al.  (1987) ,  Wyatt  and  Chang 
(1994) ,  and  Yang  et  al.  (2000b) . 


Conclusion 

Four  fundamental  analytical  representations  have  been  given  for  time-harmonic  ship  waves, 
steady  ship  waves,  and  wave  diffraction-radiation  without  forward  speed  : 

1.  The  potential  and  velocity  representations  (14)- (18)  define  a  flow  within  a  free-surface 
potential-flow  region  in  terms  of  the  flow  at  a  boundary  surface  :  the  potential  representation  de¬ 
fines  the  velocity  potential  (j)  in  terms  of  boundary  values  of  (t>  and  its  normal  derivative  d<t>/dn ,  and 
the  velocity  representation  obtained  in  Noblesse  (2000)  defines  the  velocity  u  in  terms  of  boundary 
values  of  u.  These  boundary-integral  representations,  given  here  for  deep  water,  can  be  extended 
to  uniform  finite  water  depth.  The  velocity  representation  defines  u  within  a  flow  domain  in  terms 
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of  source  and  vortex  distributions  with  strength  equal  to  the  normal  and  tangential  components  of 
u  at  the  boundary  surface.  Thus,  the  velocity  representation  does  not  involve  the  potential  (p ,  and 
can  be  used  to  couple  a  nearfield  flow  calculation  method  based  on  the  Euler  or  RANS  equations 
(for  which  a  velocity  potential  cannot  be  defined)  and  a  farfield  potential  flow  representation.  The 
velocity  representation  defines  u  in  terms  of  first  derivatives  of  the  Green  function  G  (whereas  u 
can  only  be  obtained  from  the  potential  representation  via  analytical  differentiation  of  p ,  which 
involves  second-order  derivatives  of  G ,  or  via  numerical  differentiation  oi  (j)) . 

2.  A  practical  representation  of  the  super  Green  function  Q  defined  by  (23)  has  been  given.  The 
dispersion  functions  D,  D'  and  the  amplitude  function  A  in  (23)  are  generic.  Thus,  the  super  Green 
function  (23)  is  representative  of  a  broad  class  of  dispersive  waves  (including  water  waves  in  finite 
uniform  water  depth  and  internal  waves  in  a  density-stratified  fluid)  generated  by  an  arbitrary 
distribution  of  singularities  (e.g.  source,  dipole,  vortex  sheets) .  Expression  (31)  provides  a  formal 
decomposition  of  the  singular  double  Fourier  integral  (23)  in  terms  of  a  wave  component  and  a  local 
component.  The  wave  component  in  this  decomposition  is  defined  by  (26)  in  terms  of  single 
Fourier  integrals  along  the  dispersion  curves  associated  with  the  dispersion  relation  D  =  0 .  The 
local  component  is  given  by  (29)  and  (30) .  The  single  Fourier  integrals  (26a)  and  (30a)  and  the 
double  Fourier  integral  (29b)  can  be  evaluated  in  a  straightforward  and  very  efficient  manner  if  the 
amplitude  function  A  is  not  rapidly  oscillatory  (as  is  the  case  in  practice  if  the  boundary  surface  is 
divided  into  patches) .  The  single  Fourier  integrals  (26a)  and  (30a)  can  obviously  be  combined. 

3.  The  Fourier-Kochin  approach  and  the  Fourier  representation  of  super  Green  functions  have 

been  used  to  express  the  potential  (p  and  velocity  u  defined  by  the  boundary-integral  representations 
(14)-(18)  as  the  sum  of  a  wave  component  and  a  local  component  that  represents  a 

nonoscillatory  local  flow  disturbance.  The  wave  component  is  given  by  the  single  Fourier  integrals 
(33) ,  where  summation  is  performed  over  all  the  patches  that  represent  the  surface  EUFUS^,  and 
the  wave-spectrum  function  is  defined  by  (21)  and  (22)  in  terms  of  distributions  of  elementary 
waves  over  patch  p.  Thus,  the  farfield  and  nearfield  waves  generated  by  a  velocity  distribution  at 
a  boundary  surface  are  defined  in  terms  of  remarkably  simple  Fourier  superpositions  of  elementary 
waves  over  the  boundary  surface.  Specifically,  time-harmonic  ship  waves  are  defined  by  (34)  with 
(36)  and  (37) ;  time-harmonic  waves  in  the  special  case  F  =  0  (no  forward  speed)  and  steady  ship 
waves  (zero-frequency  limit  /  =  0)  are  respectively  given  by  (38)  and  (39) .  The  Fourier-Kochin 
wave  representation  has  been  obtained  from  the  potential  representation  in  Noblesse  and  Yang 
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(1995)  and  from  the  velocity  representation  in  Noblesse  (2000) . 

4.  The  Fourier-Kochin  approach,  the  Fourier  representation  of  super  Green  functions,  and  the 
Rankine-Fourier  decomposition  process,  expounded  here  for  free-surface  Green  functions,  show 
that  the  potential  (p  and  velocity  u  defined  by  the  boundary-integral  representations  (14)-(18) 
can  be  expressed  as  in  (48) .  The  nearfield  flow  representation  (48)  extends  the  Fourier-Kochin 
representation  of  waves  given  by  (33)  and  (21)-(22) .  Indeed,  the  wave  component  (p^,  in 
(48a)  corresponds  to  the  nearfield  waves  given  by  the  Fourier-Kochin  wave  representation.  The 
other  components  in  (48a)  represent  local  flow  disturbances.  These  local-flow  components  consist 
of  Rankine  components  and  Fourier  components  that  are  respectively  given  by  distributions  of 
elementary  Rankine  singularities  and  elementary  waves  over  the  boundary  surface  E ,  the  boundary 
curve  r,  and  (eventually)  the  mean  free-surface  plane  (if  PT^O) .  The  Fourier  components  are 
given  by  single  Fourier  integrals  along  the  dispersion  curves  and  double  Fourier  integrals,  as  in 
(30)  and  (29b)- (29c) .  The  integrands  of  the  double  Fourier  integrals  defining  cp^ ,  are  finite  at 
A:  =  0  and  continuous  everywhere  in  the  Fourier  plane,  and  decay  rapidly  as  k  oo.  The  Rankine 
and  Fourier-Kochin  representation  (48)  of  nearfield  flows  only  involves  elementary  functions  and  is 
remarkably  simple. 

The  foregoing  four  analytical  representations,  which  are  the  major  results  underlying  the  Fourier- 
Kochin  theory,  provide  a  new  mathematical  basis  for  computing  free-surface  flows  about  nonlifting 
and  lifting  bodies  using  free-surface  Green  functions.  Practical  applications  of  the  Fourier-Kochin 
representation  of  waves  are  given  (for  the  case  of  steady  ship  waves)  in  Guillerm  and  Alessandrini 
(1999)  and  Yang  et  al.  (2000a, b) .  Specifically,  the  Fourier-Kochin  representation  of  steady  ship 
waves  is  coupled  with  nearfield  calculations  based  on  the  RANS  and  the  Euler  equations  in  Guillerm 
and  Alessandrini  (1999)  and  Yang  et  al.  (2000a) ,  respectively,  and  is  applied  to  the  design  of  a  wave 
cancellation  multihull  ship  in  Yang  et  al.  (2000b) .  The  representation  of  super  Green  functions 
has  also  been  used  in  Noblesse  and  Chen  (1995)  to  obtain  a  simple  representation  of  the  Green 
function  for  wave  diffraction-radiation  at  low  forward  speed  in  terms  of  the  Green  function  for  wave 
diffraction-radiation  without  forward  speed.  The  Fourier-Kochin  theory  is  currently  being  applied 
to  time-harmonic  ship  waves. 

Yang  et  al.  (2000a,b)  show  that  calculations  based  on  the  Fourier-Kochin  approach  are  remark¬ 
ably  efficient.  This  efficiency  stems  from  the  fact  that  evaluation  of  the  wave  component  (p^ , 
and  the  local  components  cp^ ,  Uq  and  (p^ ,  in  the  Rankine  and  Fourier-Kochin  nearfield  flow 


34 


representation  (48a)  involves  two  0{N)  operations  :  (i)  evaluation  of  the  wave  and  local  spectrum 
functions  5^  and  S^,  which  is  a  operation  where  g^anfjg  foj-  ^he  number  of 

panels  used  to  approximate  the  boundary  surface,  and  (ii)  evaluation  of  the  Fourier  integrals  (26a) , 
(30a)  and  (29b) ,  a  0(iVP®^^^®)  operation  where  is  the  number  of  points  at  which  the  flow 

is  computed.  Thus,  evaluation  of  the  free-surface  components  in  the  Rankine  and  Fourier-Kochin 
representation  involves  0(ArP^^®^®)+0(A^P“^^®)  operations,  instead  of  0(1VP^’^®^®)  x  0(iVP“^^®) 
operations  for  the  classical  Green  function  method.  This  property  directly  follows  from  the  identity 

gA:(C  +  2)-j[a(C-a:)+/3(7;-y)]  _  ^ki; -i(a^  + Prj)  z+i{ax  + Py) 

which  is  used  in  the  Fourier-Kochin  representation  of  free-surface  effects,  but  is  ignored  in  the 
classical  Green  function  approach.  Evaluation  of  the  Rankine  components  (j)^,  and 
in  (26a)  is  a  0(ArP^^®^®)  x  0(ArP®^^*®)  operation,  which  indeed  consumes  the  largest  share  of 
computing  time  for  the  relatively  large  values  of  and  that  are  required  at  low 

speed  and/or  high  wave  frequencies  {N  is  approximately  proportional  to  1/F^  and  ) . 

Although  the  free-surface  boundary  conditions  (3)  correspond  to  Kelvin  linearization  about  the 
uniform  stream  opposing  the  ship  speed,  the  nonhomogeneous  forcing  term  p  in  (3)  can  account  for 
linearization  about  a  more  complex  base  flow,  e.g.  flow  about  a  ship  in  the  presence  of  a  rigid  free 
surface.  Thus,  the  potential  and  velocity  representations,  the  Fourier-Kochin  wave  representation, 
and  the  Rankine  and  Fourier-Kochin  nearfleld  flow  representation,  can  be  used  within  the  context 
of  various  free-surface  linearizations. 
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Fig.  1  Dispersion  curves  for  T — 0.2, 0.25  and  0.3  Fig.2  Dispersion  curves  in  the  Fourier  planes 
in  the  Fourier  plane  (a,  l3)Flf  (a, ^)//2  and  (a, 


Appendix  1 


An  elementary  verification  of  (2a)  is  given  here.  Expression  (2a)  yields 


1  1 

e-kiC-z) 

k>4 

/  d/3  da  j  < 

>£ 

if  < 

r  27r  j 

—oo  J~oo  ^ 

^  e^U-z)  J 

[c<d 

It  follows  that  (1/r)  =  0  if  C  7^  >2: .  The  flux  through  the  two  planes  C  =  ^  i  0  is  given  by 

[a(-l/r)/5C]^lo  =  -  /  d/3  =  Air  -  x)  6{t)  -  y) 

TT  J~oo  J—oo 

Thus,  (2a)  represents  the  velocity  potential  created  at  by  a  unit  source  located  at  x. 

Appendix  2 

Expressions  (11a)  and  (2b)  show  that  that  the  free-surface  component  in  (10a)  satisfies  the 
Laplace  equation.  Expressions  (10),  (2a)  with  C>z,  and  (11a)  yield 

1  roo  roo  n  c 

47rG=lim  —  /  d/S  [  da{2A^e'^^  -  ^  e-'^^ 

e-+0  27r  y_oo  7-00  k  ‘  De 

At  the  free-surface  plane  C  =  0  we  then  have 

G(^-  fG  +  i2TG^  +  F^G^^-e{ifG-\-FG^)  = 

1  /-oo  />00  ^ 

lim  —x  /  d/d  da  —  £  with 

A  =  2kA^+D,-{f-Fa){f-Fa  +  ie){2A^- DJk ) 

Expressions  (4)  and  (lib)  yield  A  ~  ieDD'/k  as  e—^O.  Thus,  K  =  0{e)  if  D^^O  and  h  =  0{e^)  if 
£)  =  0  as  e— >0. 


Appendix  3 


The  representation  (28)  holds  if 


■>w 


foo  roo  C  p  poo 

■  d(3  da  Y,EjAj-  =  J2  ds  dnEjAj 

J-OoJ-OC  py_Q  jy_QJD=oJ-00 


D 


Here,  £  =  +  ,  and  ds  and  dn  are  the  differential  elements  of  length  along  a  dispersion 

curve  D  =  0  and  its  normal  VD .  The  relation  dn  =  dD/\\VD\\  and  expression  (28b)  show  that  the 
foregoing  expression  is  identical  to  (12)  in  Noblesse  et  al.  (1999) . 
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